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We generalize a previously proposed renormalization and computation scheme 
for nonequilibrium dynamics to include finite temperature and one-loop selfcon- 
sistency as arising in the large- N limit. Since such a scheme amounts essentially 
to tadpole resummation, it also includes, at high temperature, the hard mass 
corrections proportional to T 2 . We present some numerical examples at T = 
and for finite temperature; the results reproduce the essential features of those 
of other groups. Especially we can confirm a recently discovered sum rule for the 
late time behaviour. 
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1 Introduction 

Nonequilibrium dynamics in quantum field theory has become, during the last 
years, a very active field of research in particle physics [1-5], in cosmology [6-18] 
, and in solid state physics [[H|. The application of the general formalism of 
nonequilibrium quantum field theory [pOJ ^H has been limited up to now mostly 
to spatially homogeneous systems (see however f22fl ). The typical field theoretic 
system considered in this context consists in a classical, spatially constant field 
4>(t) (Higgs, inflaton, condensate) and a quantum state of fluctuations of the same 
or another field, chosen initially as a Bogoliubov-transformed vacuum state or a 
thermal state. The classical field is started with an initial value away from a local 
or global minimum of the classical or effective action. The time development of 
this coupled system is then studied including the back reaction of the quantum 
field in one-loop, Hartree or large- N approximations. The basic equations of mo- 
tions have been derived by several groups and a considerable number of numerical 
studies has been performed. It has been found that the system is far from being 
Markovian, showing a long-time memory. The typical late-time behaviour is a 
stationary oscillation of the field <p(t) and of the quantum fluctuations [^, |24|. |3[. 



Contrary to naive expections the classical field does not come to rest and the 
quantum fluctuations do not thermalize in the various approximations that have 
been studied. 

The question of a concise and covariant renormalization may not be the most 
urgent one in some of the above-mentioned contexts. Logarithmic corrections in 
the fluctuation integral, typically with coefficients A/167T 2 , are of course small 
in most cases; covariance may not be important in solid state applications; and 
in the case of effective theories like the sigma model for disordered chiral con- 
densates renormalization may become replaced by physical cutoffs. Also, due to 
parametric resonance, suppressed by the back reaction, quantum fluctuations de- 
velop most strongly in the small momentum region; so a cutoff chosen with taste 
will do it for practical purposes. Nevertheless, in an expanding field of research 
as nonequilibrium quantum field theory one should make sure that things can 
be done properly. In a - still unexplored - GUT phase transition the coupling 
will not be as small as for the inflationary models being investigated at present. 
Renormalization of the large- N approximation certainly will be important as a 
basic step if one tries to include a real rescattering of quantum fluctuations; such 
a rescattering - not only through the back reaction with the classical field - is 
presumably an important ingredient for understanding thermalization. A subject 
that arises in close connection with renormalization are singularities in the time 
variable due to initial conditions; it is a rather fundamental subject for nonequi- 
librium dynamics. We have addressed this problem recently |25| . Finally, from 



a practical point of view, our renormalization procedure also implies improving 
considerably the convergence of momentum integrals, and the gain in computing 
time can be of importance when studying more complex systems. 
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The basic method for the perturbative expansion which we will use here has 



been developed in f27j. It is based essentially on a standard resolvent expansion 
of the progagator. Most other groups have used the eikonal expansion which 
seems a natural choice in the presence of an oscillating background field and has 
its merit in allowing the study of adiabatic properties ||. However, renormaliza- 
tion becomes more cumbersome as the relation to the Feynman graph expansion 
becomes more remote. Furthermore, eikonal expansions cannot be generalized to 
coupled fields (one would need time-ordered exponentials, not suited for numeri- 
cal computations). 

We will present in the next section the O(N) model, the nonequilibrium equa- 
tions of motion for the classical field and the fluctuation modes, the energy density 
and the pressure in large- iV approximation and in unrenormalized form. In order 
to prepare the discussion of renormalization we give the basic equations for a 
perturbative expansion of the mode equations in section 3. Large-N renormaliza- 
tion, including finite temperature corrections is derived in section 4. Section 5 is 
devoted to the high temperature limit of the model. Some numerical experiments 
are presented in section 6. A resume and conclusions are given in section 7. 

2 Formulation of the large- N equations 

We consider the O(N) with the Lagrangian f] 

C = l -d^d^ - ^mW - ^(0> 4 ) 2 (2-1) 

where (p' l ,i = 1,..,N are N real scalar fields. The nonequilibrium state of the 
system is characterized by a classical expectation value which we take in the 
direction of <p^. We split the field into its expectation value <j) and the quantum 
fluctuations ip via 

<^(x, t) = S^VN^t) + ^(x, t) . (2.2) 

In the large N limit one neglects, in the Lagrangian, all terms which are not of 
order N. In particular terms containing the fluctuation ip^ of the component 
4>n are at most of order yN and are dropped, therefore. This is in contrast 
to the Hartree approximation where the fluctuations of <p^ are included. The 
fluctuations of the other components are identical, their summation produces 
factors N — 1 = ]V(1 + 0(1/N)). Identifying all the fields ?/>i, ..ipN-i as ip the 
leading order term in the Lagrangian then takes the form 

C = N(C (f> + C^ + C I ) , (2.3) 



We deviate from the usual convention of introducing the interaction term with a factor 

| and, even worse ^j- 



A/8iV in order to avoid a plethora of fractions -| and, even worse "t , m the subsequent 



formulae. 



with 






\d,<pd^ - ~mV - ^0 4 , 



(2.4) 
(2.5) 
(2.6) 



where ip 2 is to be identified with J2 ip l ip % /N . 

We decompose the fluctuating field into momentum eigenf unctions via 



V>(x,t) 



d 3 fc 



(2vr) 3 2^ 



a k C/ fc (t)e ikx + 4[/ fc *(t)e- 



ikx 



(2.7) 



with u>ko = y^o + ^ 2 - The mass m will be specified below. This field decompo- 
sition defines a vacuum state as being annihiliated by the operators a^. 

The equations of motion for the field <f>(t) and of the fluctuations Uk(t) have 
been derived in this formalism by various authors || |26|, §. In addition to the 
large-iV Lagrangian, Eqs. 
fluctuations, rules like 



one uses, on averaging over the quantum 



m 2 

d{ip 2 ) 2 

dip 

d 2 (ip 2 ) 2 
dip 2 



m 2 , 

AiP{ip 2 ) 
A(iP 2 ) 



or 



(2.8) 
(2.9) 

(2.10) 



which follow at large N from the identification ip 2 ~ J2 ip l ip l /N. 

We include in the following the counter terms that we will need later in order 
to write the renormalized equations. So the equation of motion for the field <p 
becomes 



4>(t) + (m 2 + 5m 2 )<P{t) + (A + 6\)<f)(t) <p 2 {t) + F(t, T) 







(2.11) 



Here T{t, T) is the divergent fluctuation integral; it is given by the average of 
the fluctuation fields defined by the initial density matrix. For a thermal initial 
state of quanta with energy u k0 = 



k 2 + m 2 ) it is given by 

d 3 k ., f3u k0 



*«.1 = <*V OH/a-s-coth 



(27T) 3 2cu fc0 
The mode functions satisfy the equation: 



U k (t)\' 



£ + <4W 



U k {t) = 



(2.12) 



(2.13) 



and the initial conditions 

t4(0) = 1 ; Lr fc (0) = -iu k0 . 
The time dependent frequency u>k(t) is given by 

ul = k? + M 2 {t) 
with the time dependent mass 

M 2 (t)=m 2 + (X + 6X)[(f) 2 (t)+F(tj . 
As in our previous work we rewrite the mode equation in the form 

d 2 



dt 2 



+ u 2 k0 



U k (t) = -V(t)U k (t) 



(2.14) 
(2.15) 
(2.16) 

(2.17) 



whereby we have defined the time-dependent potential V(t) = M. 2 {t) — .M 2 (0); 
we further identify m = .M(O) as the "initial mass". The classical equation of 
motion also can be rewritten as 



0(£)+A4 2 (£)0(£) = O 



(2.18) 



which is of the same form as Eq. ( |2.13| ) with k — 0, the classical field also is 
referred to as "zero mode" in [ij. The average of energy with respect to the initial 
density matrix is given by Q 



£ 



-0 2 (i) + -(m 2 + <5m 2 )0 2 W + 



X + 5X 



4 (t) + SA 



+ 



d 3 k 



^^coth^jil^^I^I^WI 2 



(2.19) 



X + 5X _,, 



4 



T 2 {t : T) . 



It is easy to check, using the equations of motion ( |2.18| ) and ( p,13|) , that the energy 
is conserved. The energy density is the 00 component of the energy-momentum 
tensor. The average of the energy momentum tensor for our system is diagonal, 
its space-space components define the pressure which is given by 

d 2 



P 



<t> 2 {t)-E + A 



+ 



d 3 k 



(2tt) 3 2^ 



dt 2 



coth 



-TT- w fc0 + " 



(2.20) 



\Uk{t)f 



The term proportional to A is the space-space component of the "improvement" 
term A(g^ v d 2 — d^d^cf) 2 for the energy momentum tensor as introduced by [^8[ 
It serves as a renormalization counter term, here. 



5 Note that twice the last term, with positive sign, is included in the fluctuation energy, since 
]{t) contains f(t,T). 



3 Perturbative expansion 



In order to prepare the renormalized version of the equations we introduce a 
suitable expansion of the mode functions. We have used this method exhaustively 
in our previous publications for the infiaton field coupled to itself |27j and to 
gauge bosons || in Minkowski-space and for the infiaton field coupled to itself 
in a conformally flat FRW-universe ||18| . All these calculations have been done 
for T = 0. The renormalization procedure does not change for T^O. Therefore 
we give here only a brief review of the perturbative expansion. For details the 
reader is referred to our previous work. 
The mode functions can be written as 



cP_ 



+ ^fe0 



with 



V(t) 



U k (t) = -V(t)U k (t) 



M 2 (t)-M 2 (0), 



k 2 + M 2 (0) 



1/2 



(3.1) 



(3.2) 
(3.3) 



(for the definition of M. 2 (t) see Eq. fl2.16|) ). The mode functions satisfy the equiv- 
alent integral equation 



U k (t) = e-^ ot + J dt'A k , Tet (t - t')V(t')U k (f) . 



with 



Afc, re t(^ — t ) 



1 



U k Q 



6{t-t') sin (ukoit-t')) 



For U k (t) we choose the following ansatz 

U k (t) = e -^\l + f k (t)) 



(3.4) 



(3.5) 



(3.6) 



to separate U k {t) into the trivial part corresponding to the case V(t) = and a 
function f k (t) which represents the reaction to the potential. f k (t) satisfies the 
differential equation 

f k (t) - 2«W fc (<) = -V(t)(l + f k (t)) , (3.7) 

with the initial conditions /fc(0) = f k (0) = or the equivalent integral equation 



f k (t) = Jdt'A kjiet (t - t')V(t')[l + f k (lf)]e^(t~t') 

o 



(31 



We expand now /&(£) with respect to orders in V(t) by writing 

/*(*) = fi 1) (t)+fiy)+fi 3 \t)+--- (3.9) 

= fi 1] (t)+fi 2 \t), (3.10) 

where f k n (t) is of nth order in V(t) and f k (t) is the sum over all orders beginning 
with the nth one: 

oo 

fi n \t) = Y,ft\t). (3.11) 

l=n 

The function f k (t) is identical to the function fk{t) itself which is obtained by 

solving ( |3.7| ). The function f k (t) can be computed by using the differential 
equation, via 

/f(t) - 2icWf (*) = -V(t)/F(t) , (3.12) 

or by iteration via 

t 

ff\t) = fdt'A kjIet (t - OVCO/^COe** ^ • (3-13) 



This iteration has the advantage for the numerical computation that it avoids 

computing f k via the small difference fj. —fl ■ However, the integral equations 
are used as well in order to derive the asymptotic behaviour as uj^q — > oo and 
to separate divergent and finite contributions. The leading orders of /&(£) are 
discussed in detail in |18, ^7], ^| at full length and we do not want to repeat it 



here. In this work we are more interested in the effects of the finite contributions 
at finite temperature and in the self consistent solving of the large- N limit. 

4 Renormalization 

We use the expansion and the definition introduced in the previous section in 
order to in order to single out the divergent terms from the fluctuation integral; 
we have 



M 2 (t) = m 2 + <5m 2 + (A + 5A){0 2 (t) + /_ 1 (m o ,T) (4.1) 

-J_ 3 (mo,T) [M\t) -M 2 (0)} + .F fin (£,T)} , 
where the finite part of JF(t, T) can be written as 

^■ t) = /p|k4/ di ' cos(2 ^ (i - t ' ,)v(i ' )coth ^ 2 
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and where the divergent integrals are defined as 

r d 3 k / 2 \ 

/d?k / 2 \ 

(2^r (l + 55^1) = M™o) + S-s(mo, T) .(4.4) 

The integrals i_fc(m ) are those which occur in the renormalization at T = 0. 
Their dimensionally regularized form will be given below. The additional tem- 
perature dependent terms S_fc(mo,T) are finite. They are defined as 

S - 3(m °' T) = 7 (2^)3 2^(^0-1) - ^ 

We derive some useful explicit expressions for these integrals in Appendix A. 

It is convenient to include these finite terms into the definition of J 7 & n (t,T). 
Then the time dependent mass takes the form 

M 2 {t) = m 2 + 6m 2 + (A + 6X) [<f(£) + J_i(m ) - I- 3 (m )V(t) + P^(t,T) 

(4.7) 
with 

F &n (t,T) =S_ 1 (mo ) T)-V(i)S_ 3 (mo,T)H-JF fln (i,T) . (4.8) 

The time dependent mass (|4.7| ) contains both renormalization constants 8m and 
5X. Furthermore, its definition by this equation is implicit, A4 2 (t) appears also 
on the right hand side of Eq. (|4.7|) in V(t). 

We now set out to fix the renormalization counter terms in such a way that the 
relation between the time- dependent mass and 0(t) becomes finite. An additional 
constraint derives from the requirement that the renormalization counter terms 
should not depend on the initial condition but only on the parameters appearing 
in the Lagrangian, i.e., A and m. For the simpler case of the one-loop equations 
this has been achieved |27| . 

We first determine 5\ by considering the difference 

V(t) = M 2 {t)-M 2 {0) (4.9) 

= (A + 8X) {<f) 2 (t) - 2 (O) - I_ 3 (m )V(t) + -F fin (t, T) - ^ fin (0, T)} 

or 

V(t)[l + (A + £\)i_ 3 (ro)] = (A + 5A){0 2 (t)-0 2 (O)-[/_ 3 (mo)-/- 3 (m)]V(t) 

+^ fln (t,T)-^ fin (0,T)} . (4.10) 



We now 


require 




X + 5X 






1 + (A 


+ <JA)/_ 


j(m) 


Solving 


with respect to 5X 


we fine 


[ 








5X-- 


A 2 /-; 
1-XI 


,(m) 

-3("l) 


Insert in 


g this relation into 


((4.1(1) 


we have 





A. 



(4.11) 



(4.12) 



V(t) = A {0 2 (t) - 2 (O) - [J_ 3 (m ) - I- 3 (m)} V(t) + F &n (t, T) - ^(0, T)} 

(4.13) 



or 



V(t) 



A 



1 + A [/_ 3 (mo) - /- 3 (m)] 



2 (t) - 2 (O) + ^ fln (t,T) - ^ fin (0,T) 



(4.14) 



This is a finite relation for the potential V(t) since the difference [i_ 3 (m ) 
I_3(m)] is finite. Using dimensional regularization 



I-3(m ) 

and therefore 



d 3 A; 



(2vr) 3 4u;^ 



reg 



16vr 2 



'2 , 47T/X 2 

- + ln— V 



m 



7 



(4.15) 



7_3 (mo) - 1-3 M 



1 jm 2 

16it 2 I mo 



(4.16) 



We now go back to equation ( |4.1| ) which we take at the initial time t = 0: 



m r 



M 2 (0) = m 2 + 5m 2 + (X + SX) 2 (O) + I_i(m ) + :F fin (0,T) . (4.17) 



This is an implicit relation between m and 0(0) which, however, contains still 
the infinite quantities SX,Sm and i_i(m ). In order to proceed we note the fol- 
lowing explicit relation between I_\ and I_ 3 which follows from the dimensionally 
regularized expressions for these quantities 



d 3 k 



m: 



o 



(2vr) 3 2^ 



16vr 2 

,2 



+ In — y- -7 + 1 



-m /_ 3 (m ) 



Therefore we can rewrite ( 4.17| ) as 



rrin 



m 2 + bra 2 + (A + 5X) 



2 (O) - m 2 /_ 3 (m ) 



m? 



16tt 2 



16vr 2 



(4.18) 



^fin(0,T) 



(4.19) 



or 



m 2 [1 + (A + 5A)i_ 3 (m)] = m 2 + 5m 2 



(4.20) 



+ (\ + 5\) 



rrin 



2 (O) - m 2 (I. 3 (m ) - 7_s(m)) - -^- + ^ fin (0, T) 



16tt 2 



We now require the factors A + 5\ and [1 + (A + 5A)i_ 3 ] to cancel on account of 
( }4.11|) , so as to obtain a relation between finite quantities. This is obviously the 
case if 

m 2 + bra 2 = m 2 [1 + (A + 5A)i_ 3 (m)] + (A + 8\)p 

since then 

^2 



///n 



m 2 -m 2 = A 2 (O) - m 2 [/_ 3 (m ) - I. 3 (m)] - — \ + p + ^(0, T) 



167T 2 



(4.21) 
(4.22) 



( |4.21D fixes bra 2 as 

5m 2 = (A + 5 A) 
Setting p = we have 



m 2 I- 3 {m) + p 



A 



m 2 I^ 3 (m) + p 
1 - A/_ 3 (m) 



5m 



Am 2 I_3(m) 
1 - A/_ 3 fm) 



(4.23) 



(4.24) 



and 

m 2 - ra 2 = A J0 2 (O) - m 2 [/_ 3 (m ) - I_ 3 (ro)] 



Til ~ 1 

T^ + ^(0,T)|. (4.25) 



This corresponds to the MS subtraction. Another natural choice is p = m 2 /167r 2 ; 

then 

m 2 /_ 3 (m) + m 2 /167r 2 AJ_i(m) 

(4.26) 



Sra 2 = A- 



AJ_,m 



1 - AI_,(m 



and 



rn 



l-ra 2 = A 2 (O) - m 2 [I_ 3 (m ) - /_ 3 (m)] 



2 2 

m,Q — ra 



16vr 2 



+ ^ fin (0,T) . (4.27) 



This choice is analogous to the one in [27| for the one-loop equations and will be 
used in the following. 

The "gap equation" (|4.22| ) and the renormalized definition of the potential 
( 4.14 ) constitute, along with the equations of motion the basic renormalized 
equations for the self consistent large- N dynamics. 

The gap equation has to be solved at t = and determines the relation 
between ra and 0(0). For later times we have 



M 2 (t) = m 2 + V(t) 

= ml + CX U 2 (t) - 2 (O) + J- fin (t, T) - ,F fin (0, T) 



(4.28) 



with 



C 



1 + 



16vr 



In 



m 



rrin 



(4.29) 



Since the gap equation can be cast into different forms we can obtain several 
equivalent forms of this equation. Solving the gap equation for 2 (O) we find 



AC0 2 (O) = -ml + C 



m 



A 



so that 



M\t) = C 



m 2 + A (f) 2 {t) 



rrin 



16tt 2 



rrin 



16vr 2 



p-Jfi n (0,T) 



P + Ffin(t,T) 



(4.30) 



(4.31) 



Having obtained a finite relation between <fr(t) and Ai(t) the equations of motion 
for the classical field (f)(t) and for the modes Uk(t) are well-defined and finite. 

Here we have chosen to include the corrections of leading order, proportional 
to £_i(m ,T), into the finite part of the fluctuation integral. These terms are 
important at high temperature; they appear in the gap equation via ^fi n (0, T) = 
£_i(m ,T) ~ T 2 /12. Omitting some terms of order \/16tt 2 the gap equation 
l.22|) becomes 

, (1 _m 2 + A0(O) 2 + 4r 2 . (4.32) 



m 2 , ~ 



12" 



Therefore, at high temperature the mass circulating in the loop is dominated by 
the "hard" AT 2 term. 

We will need in the following the fluctuation integral J-(t, T) which is and will 
remain divergent. We need an expression in which these divergencies appear in 
explicit form. We use 



J r (*,T)=I_ 1 (mo)-i-3(m ) M 2 (t) - M 2 (0) +^ fin (i,T) 



(4.33) 



and insert the expression for A4 2 (t) we have just derived. Using the gap equation 
and some reshuffling of terms we obtain 



F{t,T) 



rrin 



CAJ_ 3 (m o )0 2 (t) - m 2 CI. 3 (m ) 



16vr 2 
+ C T7T^ J -3( m o)K - m 2 ) + C{1 - A/_ 3 (m)) P 6n {t, T) . 



(4.34) 



5 Renormalization of energy and pressure 

The expressions for the energy density and for the pressure have been given in 
section 2. Apart from the renormalization counter terms which we have already 
fixed in renormalizing the equation of motion, two new counter terms appear, the 
"cosmological constant" term 8h in the energy density and the "improvement 
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term" Ad 2 {(f) 2 + (ip 2 ))/dt 2 in the pressure. These terms must suffice for rendering 
the expressions for energy density and pressure finite. 



\^ 2 (t) + imvw + \m 



We start with the expression ( 2.19 ) for the energy which we rewrite as 
£ = 

with 



^^(t,r) + >V(^ ^oVu ''A, (T).i) 



i 

— < 

2 



T 



_H_ cotll /^(i l ^ (i)P+ ^ w|c/tW|2 



In the latter expression we split the Bose factor as before 

f3u) k0 . 2 



coth : 



1 + 



2 ' e /3"fco _ i • 

The integrations involving the second term are finite, we define 



A£ fl (£,T) 



(X eM!-i &^| 2 + ^l^| 2 



(5.2) 



(5.3) 



(5.4) 



Those involving the first term have been discussed in P7[ . Following this discus- 
sion we can decompose the integral via 



£ & {t,0) = h(m ) + iv(t)/_i(m ) - ~V 2 (t)J_ 3 (m ) +£i4i a (*,0) 



(5.5) 



withg 



£fl,fin(£, 0) — - 



rf 3 A; 



1. : 



2 7 (27r) 3 2cj fc0 12 



ol/j 



(1)|2 



+ 



V(t) 



r(l) 



2Re/r J + |/i 



(1)|2 



+ 



V 2 (t) 



So; 



A0 



(5.6) 

We denote the sum of both finite contributions as £s.fi n (t, T). The expression for 
the energy now takes the form 

£ = l0 2 + I m V + ^ 4 + ^, fi n(t,T) + / 1 (mo) + ^V(t)/_ 1 (mo)-^V 2 (t)/_ 3 (mo) 



■^(^ + i*my + ^ + W. 



(5.7) 



In addition to the divergent integrals I n (mo) and the counter terms further diver- 
gencies are contained in the fluctuation integral JF(t, T); these are given explicitly 



6 We have overlooked in 27 that in the expression for the fluctuation energy given there 
two terms cancel on account of a Wronskian identity, given in Q, Eq. (73). They did so, of 
course, in the numerical calculations. 



11 



in Eq. ( |4.34| ). The analysis of Eq. (|5.7|) , after inserting the explicit expressions 
for V(t) and J-(t, T), becomes rather cumbersome. To give an outline of the typi- 
cal algebraic manipulations we consider explicitly the coeffi cents of 4 . Collecting 
everything except the bare \<j) 4 term we find that 4> (t) is multiplied by a sum of 

divergent terms 

5X 

4 

We use the abbreviations xq = 7_3(mo) and x = 7_3(m), so that 



±±^C 2 X 2 x 2 - -C 2 x X 2 
4 4 



(51 



C 



1 



X(xq — x) 



(5.9) 



and (see Eq. ( gig )) 



sx 



X 2 



1 — Ax 
One finds that all divergent quantities combine into the finite expression 



(5.10) 



-X 2 C(x — Xq) 



X - £ J% 



64tt : 



m 



AA 



(5.11) 



so the correction to the 4 term in the energy becomes 

AA .„ 



(5.12) 



Collecting similarly all terms proportional to <f) 2 one finds that the correction to 
the mass term becomes finite as well, explicitly 



■Am 



2/2 



A 



32tt 2 



-C 



m 



rrin 



m 2 In 



m 



(5.13) 



There are further time-dependent terms proportional to Jg n (t, T) and jF fin (t, T) 
and constant terms. The divergent parts of the latter ones can be absorbed into 
SA, the coefficient of the term linear in T^ a vanishes and the quadratic one has a 
finite coefficient. The counter term 5A can be chosen independent of mo: 



5A 



m 



4(1 - Ax) 



x + 



A 



57T" 



256vr^ 



(5.14) 



there remains a finite constant 



AA= -C 

4 



(x — x)m 4 H -(m 2 



i-K^ 



m 2 Q ) + 



32tt 



4 i 

r m + 



256vr 



2 2\2 

4 v m - m ) 



(5.15) 
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So the expression for the energy can really be rendered finite with counter terms 
independent of the initial condition. Explicitly we find 

A + AA, 4 



£ = -0 2 + -(m 2 + Am 2 U 2 4 

2 2 4 

+£ fl , f (t,T)-^ 2 n (t,T) + AA 



(5.16) 



We finally have to give a finite expression for the pressure, using our last free 
counter term. We write the pressure in the form 



p = 2 (t) - £ + pa (t, T) + A^ [0 2 (t) + F(t, T) 



(5.17) 



Here we have anticipated a special form of the counter term, indeed for the 
expression in brackets one can choose a priori an arbitrary Lorentz scalar, the 
additional piece of the energy momentum tensor being trivially conserved on 
account of its tensor structure d^d u — g^d 2 . Of course it has to be suited for the 
renormalization procedure. The fluctuation part of the pressure consists again 
of three parts, a divergent one, a finite one independent of the temperature and 
a finite integral involving the thermal distribution function l/(exp(u)ko/T) — 1). 



The analysis for T = has been performed in [27|. Following the discussion there 
we can write p& as 



d 3 k 



Ap t s.(t,T) is given by 

Ap )fl (t,T) = 
the T = finite part by 

Pfl,fin(£, 0) = 



m 



m 



96tt 2 48tt 2 



6 



48^ 



(2vr) 3 2u; fc0 e/*"*> - 1 



V(t) - -(/_ 3 (mo) + T^Mt) ■ 

(5.18) 

(5.19) 



CO 



w+ 3 



\Uk(t)f 




+ 



6^o 



,12^o 



2*4>, 



+ 



m? 



2M 0/ 



2Re/f(t) + |jf(t)| 2 

fdt cos2uj k (t-t') V (*') 
o 

) cos(2^ t)V(0) 



+ |/i 1) (t)| 2 - 2Re i^oPHt) +iu ko fl 1) (t)fi 1) *(t) 



(5.20) 
(5.21) 



We call the sum of both finite fluctuation integrals pafi n (t,T). Now we have to 
consider the divergent terms. We observe that V(t) is given by 



V(t) = AC^[0 2 (t)+J- fin (t,T) 



(5.22) 
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On the other hand, using Eq. (|4.34|) we have 



nr^ T) 



dt 2 



and therefore 

d 2 



-\Cx 4> 2 (t) + C(l - \x)F &n {t, T) 



d 2 



A ^ {(f)2{t) + ^ (t ' T)) = A d^ c{1 ~ Xx) r (t) + ^ fin(t ' T) . 



(5.23) 



(5.24) 



As apparent from Eq. ( |5.22|) this matches in form with the divergent term 
J_3(wio)V(t)/6. Insisting again in choosing the counter term independent of the 
initial condition we fix 



.4 



Ax 



A/_ 3 (m) 



6(1 -Ax) 6(l-AJ_ 3 (m)) 



and retain a finite term 



96tt 2 



in ™ 

TO . 



v(0- 



(5.25) 



(5.26) 



The final result for the pressure reads 



p = 2 (£)-£+p fl , fin (£,T) 



m 



r?i, 



96tt 2 48tt 2 



^V(t) 



96tt 



In ™ 
Wo, 



V(t) . (5.27) 



A further quantity of interest is the particle number density. It does not need 
to be renormalized after subtraction of the initial particle number density and is 
given by 



n{t) - n(0) 



d 3 k (3u k0 j 1 

(27F coth Ml 



\u k (t)\ 2 + ^rlUkW 



U! 



A-0 



2 



d 3 k coth Puko\f (1 Kt)\ 2 



(2tt 



4^lo 



(5.2* 



6 Numerical results 

We have implemented numerically the renormalized formalism derived in the 
previous sections. The results show essentially the same features as those found 



by other groups [30 



We have chosen several parameter sets which are displayed in Table 1. All 
parameter sets start with an initial value of the mass mo = 3. The gap equation is 
then solved for 0(0). We have chosen two parameter sets at T = 0, one with A = 1 
and one with A = 5. Another two parameter sets have finite temperature, T = 3 
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and T = 10, respectively, and A = 1. Parameter sets with smaller temperatures 
showed very little deviations from the T = case and are not presented. 

For the first parameter set the numerical results are diplayed in Figs. 1 
a-e. The classical field (pit) is seen to oscillate with an amplitude decreasing 
slowly to an asymptotic value. The potential V(t), the difference between A4 2 (t) 
and m,Q = .M 2 (0) reaches an asymptotic average of —4.0 with small oscillations. 
Classical and fluctuation energy are shown in Fig. lc. We denote as fluctuation 
energy the quantity £ fljfin (t,T) = £ flifin (£,0) + A£ fl]fin (£,T), see Eqs. fl5.6|,|5.4j) 



The remaining parts of the total energy (|5.16|) are considered as the classical 
energy. The separation is somewhat arbitrary as finite parts of the leading order 
fluctuation energy are contained in AA, Am 2 and AA. So the fact that the 
"classical energy" becomes even negative is deceptive, the classical amplitude 4>(t) 
has not decreased to zero. It has decreased roughly by a factor of 3, implying a 
decrease of energy by a factor 9. Nevertheless the production of fluctuation energy 
is important, as also seen from the particle number displayed in Fig. le. The 
energy is seen to be conserved, apart from some numerical noise in the early stage 
of evolution, due to badly convergent integrals. The pressure is seen to approach 
an asymptotic value of ~ 5, somewhat smaller than the ultrarelativistic limit of 
E/3 = 6.7. 

For the other parameter sets we display just the classical amplitude and the 
potential V(t), need to verify a sum rule (see below). For the finite temperature 
case T = 3 the particle number develops almost identically to the one for pa- 
rameter set 1, and reaches an asymptotic value of ~ 8.5. This can be compared 
with the thermal particle number density 3.29. The situation changes strongly 
in the high temperature situation, parameter set 4. There, the particle produc- 
tion, displayed in Fig. 4c, is insignificant with respect to the thermal particle 
number density which is n(0) ~ 122. The other Figures for this parameter set 
show clearly that the interaction of the classical amplitude with fluctuations is 
suppressed. 

An interesting topic which has emerged recently |J(J is a sum rule for the 
late time behaviour, which was found to be satisfied numerically with high preci- 
sion. Adapted to our notation and definitions, and generalized (naively) to finite 
temperature, it reads 

M 2 (oo) ~c(m 2 - -^K - ™ 2 ) + E_i(mo, T) ] + AC^ . (6.29) 

y 1671^ ) 2 

This value of .M(oo) is determined by the lower limit of a parametric resonance 
band, essentially it implies that the classical oscillation becomes stationary if its 
frequency, i.e. M.(t), settles in such a way as to avoid resonant excitation. We 
have verified this sum rule for our parameter sets, the left and right hand sides 
of the sum rule are compared in Table 1. For T = the agreement is excellent, 



for T^0,a case not considered in [30|, the deviations are of the order of 10%. 
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7 Conclusions 

We have derived in this paper the renormalized equations of motion, the energy 
and the pressure for the nonequilibrium evolution of a scalar O(N) model. The 
regularization has been done in a covariant way, using dimensional regularization. 
As in the case of the one-loop equations studied previously [^7| it was possible 



to fix all the renormalization counter terms independent of the initial conditions, 
though the divergent integrals appearing in the unrenormalized expressions do 
depend on the "initial mass" mo instead of the renormalized mass. As a renor- 
malization convention we have chosen a slightly modified MS scheme, it can 
of course be modified to another suitable convention like renormalization at the 
minimum of the effective potential. 

We have restricted the formalism to the case of unbroken symmetry. The 
special - and highly interesting - aspects of the broken symmetry case have been 
investigated in ||. Here it was our aim to present a general framework, which has 
to be adapted to specific physical models. In the case of the one-loop equations 
it was indeed possible to extend the formalism to nonabelian gauge theories ||. 

The formalism developed here represents at the same time a rather convenient 
computation scheme. The CPU time requirements are of the same order as the 
one for the one-loop equations. Typically, the examples we have presented took 
1-2 hours each on a small workstation. We have not attempted to meet the same 
standards in numerical precision as other groups, nevertheless our results show 
the same general features as those of other groups. This is presumably due to 
the fact that the system as such has a stable and essentially predictable late time 



behavior ||30|| , indeed our results fulfil an asymptotic sum rule formulated in this 
Reference. 



A Some thermal integrals 

In this Appendix we give, without claim of originality, some explicit expressions 
for the thermal integrals as we have used them in the numerical computations. 
In deriving these relations we have relied on the integral tables of Prudnikov, 



Brychkov and Marichev [29 



The finite temperature part of the tadpole graph, which constitutes a correc- 
tion to the mass, is given by the integral 

d 3 k 1 




(2tt) 3 c^ (e^o _ l 

E 



(A.l) 



3=0 



4j!(j + 1)! 



21n%-V(j + l)-^(j + 2) 




where T nm stands for nmo/T '. For large values of T nm (this means for small T) the 
integrand is dominated by momenta of order k ~ T. Therefore one can expand 
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Uko w.r.t. powers of m/k; the integral is then well approximated by 

E^mo, T) c* || £ yfe^Tni {l + ^ - ^T nm 2 + ^T_ 3 + (T nm 4 

(A.2) 
For T ^> mo we find directly from Eq. ( |A. 1| ) the well-known approximation 

S_ 1 (m ,T)^^C(2)T 2 = ^. (A.3) 

It yields the hard thermal loop corrections to the mass. 

The finite temperature part of the fish graph, which can be considered as a 
finite correction to the coupling constant, is given by 



T 

■*■ 1) .'. 



-I OO OO I 1 



'X?_V 71 



-2 

'nm i 

[ nra 



+ 2 X, 



For small T or large T nm we find the approximation 

2,_ 3 (m ,i ) ~ ~^2 2^ ~Y e lnm y- ~ Y nm ~ "128" m _ 1024 V^nmjj (A b J 

For large temperatures this integral behaves linear in T, more precisely 

£_ 3 ( TOo ,T)~-^-. (A.6) 

87r m 

The finite temperature part associated with the quartic divergence in the energy 
is given by the Planck formula 

f d 3 k Uko 

Si(»"o,r) = J {27r)3 ^ _ j 

As an approximation for large T nm or small T we find 

ml * _ T /7f -| r 27 j 705 2 2625 3 / a"I 

Ll -2^^ e m V2 nm i + y z nm + i28 nm + To2i nm+ \ nm n 

For large temperatures one obtains 



(A.* 



7T 



2 



S^mcT)-^ 4 . (A.9) 
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Table Captions 

Table 1: Parameter sets and sum rule. We display the parameters of the 4 
numerical simulations. The mass unit is set by m = 1; 0(0) follows from the gap 
equation (|4.22|) , .M 2 (oo) is the sum of Wq and of V(oo) read from the correspond- 
ing Figures. R.h.s. is the right hand side of Eq. fl6.29|) . 



Figure Captions 

Fig. la: <fi(t) for parameter set 1. 

Fig. lb: The potential V(t) for parameter set 1. 

Fig. lc: Classical energy (long-dashed line), fluctuation energy Se_^ n (t, T) (short- 
dashed line) and total energy (solid line) for parameter set 1. 

Fig. Id: The pressure p(t) for parameter set 1. 

Fig. le: Particle number n(t) — n(0) for parameter set 1. 

Fig. 2a: <p(t) for parameter set 2. 

Fig. 2b: The potential V(t) for parameter set 2. 

Fig. 3a: <p(t) for parameter set 3. 

Fig. 3b: The potential V(t) for parameter set 3. 

Fig. 4a: <p(t) for parameter set 4. 

Fig. 4b: The potential V(t) for parameter set 4. 

Fig. 4c: Particle number n(t) — n(0) for parameter set 4. 
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Figure 1b 
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Figure 1c 
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Figure 1e 
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